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The release of conduction-band electrons from a metal surface by a sub-femtosecond extreme 
ultraviolet (XUV) pulse, and their propagation through the solid, provokes a dielectric response 
in the solid that acts back on the photoelectron wave packet. We calculated the (wake) potential 
associated with this photoelectron self-interaction in terms of bulk and surface plasmon excitations 
and show that it induces a considerable, XUV-frequency-dependent temporal shift in laser-streaked 
XUV photoemission spectra, suggesting the observation of the ultrafast solid-state dielectric response 
in contemporary streaked photoemission experiments. 

PACS numbers: 78.47.J-, 42. 65. Re, 79.60.-i, 

The sudden release and subsequent motion of a photo- 
electron (PE) in and in front of a solid dielectric medium 
provokes collective electron excitations in the solid. The 
back-interaction of these excitations with the PE can be 
modeled as a dynamic, wave-like redistribution of elec- 
tronic density in the solid in terms of a complex- valued ef- 
fective electron-self-interaction (or "wake") potential ll- 
This density wake appears since the electron distri- 
bution in the solid cannot equilibrate during the motion 
of the released electron. The wake potential depends on 
the kinetic energy E of the PE. Its real part is due to vir- 
tual excitations of bulk and surface plasmons, while its 
imaginary part accounts for inelastic scattering and en- 
ergy loss. The dependence of the dynamic wake potential 
on the charge state and velocity of a classical particle has 
been studied extensively over several decades with regard 
to energy loss JS Q , electron-exchange and -correlation 
contributions 'Vj, electron emission in ion-surface colli- 
sions ^8(1, and electronic self-interaction effects on photo- 
electron spectra Q . While these examples emphasize the 
influence of the solid's dielectric response, they do not re- 
solve the ultrafast electronic response in the condensed- 
matter-plus-charged-particle system in time. Owing to 
significant progress in laser technology over the past 
decade, sub-femtosecond XUV pulses can now be gen- 
erated and synchronized with the primary IR laser pulse, 
allowing the time-resolved observation of the electronic 
dynamics in atoms [10] and solids [11]. Time-resolved ex- 
periments at the intrinsic time-scales of an active electron 



and the correlated dynamics of two electrons 



sion of (i) the PE's phase evolution during the streaked 
emission, (ii) the attenuation of the IR pulse inside the 
solid (skin effect), and (iii) electron transport effects in 
the solid. The phase of the PE is further affected by 
its plasmon-mediated self-interaction while moving inside 
and outside the solid, and we expect streaked photoe- 
mission experiments to help reveal these dynamic many- 
body effects in solids. In this Letter we investigate how 
these three factors affect the PE dynamics and extend 
previous theoretical studies on the streaked photoemis- 
sion from solid surfaces 1 to expose the effect 
of the dynamic plasmon response on time and energy re- 
solved PE spectra. 

We calculate the dynamic wake potential assuming 
that the released PE moves with a constant velocity 
Vz > Q along a classical trajectory towards and perpen- 
dicular to the surface and crosses the metal- vacuum inter- 
face (z = 0) at time t = Q, leading to the density p(r) = 
5{y\\)5{z — Vzt) (throughout this Letter we use atomic 
units, unless stated otherwise). The semi- infinite solid is 
modeled in jellium approximation [isl ] , and its excitations 
are described by the dispersion relations for bulk- and 
surface-plasmon excitations, ujf. = uj^ + 3kpk^/5 + 
and uj^ ^ uj"^ + V3kFiOsQ/V5 + PQ^ + Q'^/i U^^. For 
low momenta k = (Q,^^), these relations model single- 
plasmon modes with bulk- and surface-plasmon frequen- 
cies ujp — inn and a;^ = ujp/ v^, respectively, and decay 
into particle-hole excitation at high momenta through the 
or terms fc^/4 and Q^/d. n is the bulk-conduction-electron 



plasmons [14lll5| promise unprecedented sensitive experi- 



mental tests of collective electronic transport phenomena 
in solids and novel plasmonic devices [3 [l3 ■ 

Using attosecond time-resolved XUV-I-IR pump-probe 
technology in a proof-of-principle experiment, a relative 
temporal shift of 110 as (1 as —10^^^ s) between the pho- 
toemission of 4f core level and conduction-band electrons 
from a tungsten (110) surface was measured pjj]. Essen- 
tial for the correct reproduction of this shift within simple 
quantum mechanical models 18, 3] is the proper inclu- 



density, and kp — (Stt^u)^/^ is the Fermi velocity. /3 is 
determined so that the surface-plasmon dispersion rela- 
tion joins the particle-hole continuum at the same point 
as the bulk line [23|. The plasmon field of the solid is 
then given by the Hamiltonian Hq — ^, >q Wfefej^fek + 



^qWQaqaq, where 5q and ajj are annihilation (cre- 
ation) operators for bulk and surface plasmons, respec- 
tively [a 0, ■ The interaction between the PE and the 
jeUium sohd is given hy Hi — J drp{r) [(j)b{r) + (j)s{r)], 
where (/)fa(r) = X;Q,fc,>o ^k^k sin(A:^z)e**^ ""ii e(--z) + h.c. 
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and (jjsir) = X]q ^QaQe"'''^'e"^ '"" + h.c. are the bulk 
and surface plasmon fields. |i?kP = SnLUp/ (Vk'^uJk) and 
|AqP = TTul/ {SQujq) are the interaction strengths (with 
quantization volume V and surface 5) 

The wavefunction \'(p{t)) of the bulk and surface plas- 
mon field interacting with the classical charged particle 



for the model Hamiltonian H = Hq + Hi can be solved 
exactly 0, Q . The wake potential of the PE is then cal- 
culated as V,m{t) = \{ii)\t)\Hij{t)\ip{t)), where Hij{t) 
is the interaction-picture presentation of Hi. Following 
[3) 0, [13 , we obtain the real part of wake potential 
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where 6(z) is the unit step function and 7 the decay- 
width of the plasmon excitation. The first term in ([T]) 
includes bulk and the last two terms surface plasmon 
excitations. 

Figure [ija) shows V[^{z^Vz) for different PE veloc- 
ities Vz for aluminum with = 0.378 [24i |. assuming 
7 — Q.lujs- The static image potential is obtained in the 
limit Vz = 0. For Vz > 0, Vim oscillates near the metal 
surface with wavelength Af, = Trvz/ujp inside the solid and 
Xs — 2Trv z/i^s in the vacuum, and approaches — l/4z (23 | 
far away from the surface for all Vz, as expected. Equa- 
tion ([T]) underestimates the influence of the positively 
charged ion cores, and its bulk limit (« 5.8 eV) does 
not reproduce the Al Fermi energy (ep = 11.7 eV) [25| . 
We therefore add a step potential to obtain the effective 
dynamic image potential. 



U{z,Vz) = - 



Vn 
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z/a 



+ Vim{z,Vz) 



(2) 



and adjust the depth to Vb — 11.7 eV and interface- 
thickness parameter to a = 1.4 A, respectively, in order 
to reproduce Sp in the static limit for Vz—Q [Fig. [TJb)]. 

We model the metal surface as a 300 a.u wide slab 
and obtain its eigenvectors e„ and and wave functions 
ipn{z) by diagonalizing the time-independent Schrodinger 
equation (SE) 



Sni'niz) 



If 

'2dz 



M^)- (3) 



In typical streaking experiments, the XUV pulse inten- 
sity is sufficiently low, so that photoemission in the XUV 
pulse Ex{t) can be treated perturbatively. The release 
and propagation of the PE wavepacket 5ipn(z, t) from the 
state ij^niz.t) is then dictated by the SE [19| 



d 

l-g^6lljn{z,t;T) ^[HL{t) + 5U{z,Vz)]5^n{z,t]T) 

+ zEx{t + T)i;n{z,t;T). 



(4) 




FIG. 1: (Color online) (a)Dynamic wake potential V[^[z,Vz) 
for Al for different PE velocities v^. (b) Real part of U{z,Vz), 
adjusted to the Fermi energy of Al at Vz = 0. 

Hl = ■kl-i^ + AL{z,t)f + U{z,0) is the Hamilto- 
nian for the solid slab in the presence of the IR-laser 
pulse AL{z,t), T the delay between the XUV and IR 
pulses with r > if the XUV precedes the IR pulse, 
SU{z,Vz) = Vim{z,Vz) — Vimiz,0) represents the complcx- 
valued dynamic part of the PE self-interaction, and Vz — 
y^2{uJx — \£n\) is the PE speed. The evolution of the n- 
th initial state of the slab below the Fermi surface in the 
IR pulse is given by 



d 

l-g^tpniz,t) = HLit)^pn{z,t). 



(5) 



Since i?L(x)(^ -> ±00) = 0, ^ and Q are solved for the 
initial conditions Sipn{z,t — > — oo;r) = and ipn{z,t — > 
-00) = 'il;„{z)e-"'"^ . 

We represent the vector potential of the laser pulse as 
AL{t) = ^osin^ (7rt/TL)cos[a;L {t - tl/2)] for < i < tl 
and otherwise, with = 1-57 eV, intensity II = 
Alul/2 = 5 X lG"W/cm2, and pulse length tl = 8 fs. 
We assume an exponential damping of the IR laser field 
inside the solid AL{z,t) = AL{t) [e^/''^e(-z) -f e(z)] 
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FIG. 2: (Color online) Streaked PE spectra for an Al surface 
calculated with the (a) dynamic image potential SU {z, v^), (b) 
static image potential SU{z,Vz — 0) for hcux ~ 40 eV, skin 
depth Sl = 0, and mean free path A = 5 A. (c) Corresponding 
center-of-energy shifts SEcoe- The temporal shift between 
the traces SEcoe for dynamic and static image potentials is 
^ 100 as. 

with a skin depth Sl, and take a Gaussian XUV pulse 
with pulse length tx = 300 as. 

Assuming free-electron dispersion, E ~ fc^/2, the 
energy-differential photoemission probability P{E,t) — 



S'll>n{k,T] 



leads to the delay-dependent center 



of energy (COE) of the PE spectrum 
1 



COE 



(r) 



^<eF 



jdk k 5-ijjn{k,T) 



(6) 



where i5'0n(fc,T) is the Fourier transform of Sipnizjt 
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J dk Sipnik, t) 



the total 



oo;t) and Ptotir) = Y.e 
emission probability. 

In order to reveal the effect of the dynamic dielectric 
response on the streaked PE spectrum, we compare the 
results of two separate calculations for values of huJx be- 
tween 30 and 100 eV. First, we solve (jH) without 5U{z, Vz) 
and denote the results as "static" . Next, we add wake ef- 
fects and obtain "dynamic" results by including the real 
part of SU{z,Vz) in (HJ. In both calculations, we replace 
the imaginary part of dU{z, v^) by the phenomenological 
expression dUpf^{z,Vz) — —VzQ{~z)/{2X) '26'], where A is 
the PE mean-free path. 

We first present results for Sl = for which the IR 
field is completely screened inside the solid. Figure [2] 
compares the static and dynamic streaked PE spectra 
and their center of energies for hux — 40 eV. We define 
temporal shifts TsJ£j(^j^„) for static (dynamic) calculations 
relative to by fitting [l8l.[l9t 



E^Mdyn)^^^ = a + ML(r - Ar,,,(,,„)). (7) 

Fig. ^c) shows the relative temporal shift Ar^ake — 
ATdyn — ATsta ~ 100 as. Suggesting a noticeable - on the 
scale of the temporal resolution in measured streaked PE 
spectra - contribution of the dynamic plasmon response 
to the temporal shift. 
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FIG. 3: (Color online) (a) Ar^afce and Arsta as a function of 
the XUV frequency obtained by fitting Eq. ([7]) to the calcu- 
lated Ecoe{t), including all levels below the Fermi surface, 
for skin depth Sl = 0. (b) Comparison of Ar^nake, obtained 
from the streaked Ecoe{t) at the Fermi surface, and the 
Wigner delay Ar^;,g. {c)ATn,ake for different skin depths. 

Our results for Ardyn and Argta as a function of 
hijjx in Fig. [21^a) reveal that the dynamic wake poten- 
tial has a significant effect on the PE delay, especially 
for hiux <50 eV, where Ardyn develops a double- hump 
structure. The huJx dependence of Ar^jofce can be under- 
stood as due to scattering of the PE in the wake poten- 
tial dU{z,Vz). This interaction of the PE with dU{z,Vz) 
changes the phase of 5ipn{z,t), giving rise to a Wigner 
delay Ar^ We determine ArJ^^^^^^ for ^z, = by 
relating the PE position (z) = J dkz\5il)n{zTt)\'^ and ve- 
locity {v) — J dk\kStpn(k,t)\'^ at a time t ^ tl according 
to (z) = {v){t — At^). In support of this "scattering in- 
terpretation" , Fig. [3Jb) shows excellent agreement of the 
streaking delay Ar^ake with Ar^j.^ = At^^ - Ar]f„ 
for photoemission from the Fermi level. We find equally 
good agreement for emission from initial conduction band 
states below the Fermi level. For this comparison, we 
assumed that the IR field does not penetrate the solid 
(Sl = 0). 

Since the effect of the actual IR skin depth on the 
streaked s pec trum from surfaces is currently debated 

ii|,[ii,[2jt2i], we found it compelling to study the infiu- 

ence of 6l on Ar^ake- Our numerical results in Fig.[3fc) 
show that Ar^ake is strongly affected by changes in the 
IR skin depth for 6l ^ 2A, due to AC Stark polariza- 
tion of the initial states 4'n{z, t) and the PE wave packet 
Stpnizji), as well as the combined action of U{z, Vz) and 
AL{z,t) on S4'n{z,t) [l^. As contributions to the emit- 
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FIG. 4: (Color online) Relative delay Arwake as a function of 
the XUV frequency at different (a) PE mean-free paths and 
(b) surface plasmon frequencies. 

ted photocurrent are limited to PEs that are released 
within a few mean-free paths A from the surface, AT^ake 
becomes less IR skin-depth dependent for 6l ^ 2A and 
stabilizes in the limit — >■ oo. 

Figure |4] shows the dependence of Arwake on the PE 
mean free path and the surface plasmon frequency for 
Sl = 0. Increasing A by a factor of two significantly in- 
creases ATwake for fujjx ^ 50 cV, but has little influence 
at larger hujx (Fig. Hfa)). Our result that, in general, 
AT^afee(2A) ^ Ar„,n,fcp.(A), is incompatible with the in- 
terpretation [llj, 122] of the delay between photoemission 
from core and conduction-band levels in tungsten being 
due solely to the PE's average travel time in the solid 
(w X/{v)). Decreasing utg shifts the double-hump struc- 
ture to the lower hux , and thus to lower kinetic energies 
of the PEs, as expected in view of the decreased thresh- 
olds for surface and bulk plasmon excitation (Fig. |4ljb)). 

In conclusion, we have calculated the wake potential 
induced in an aluminum surface by a PE that is re- 
leased in the electric field of an attosecond XUV pulse. 
By comparing centers of energies and photoemission de- 
lays in IR-streaked PE spectra including the dynamic 
wake potential with calculations performed in the static 
limit, we find a significant contribution to the tempo- 
ral shift ATmake in photoemission from the metal con- 
duction band. This shift is due to the excitation of the 
bulk and surface plasmons in the metal during photoe- 
mission and is found to sensitively depend on the XUV 
frequency as well as on solid state characteristics, such as 
the bulk (surface) plasmon frequency, IR skin depth, and 
PE transport in the solid. The measurement of streaked 
electron spectra from dielectric solids, may thus be ap- 
plied to probe solid state characteristic, in particular, the 
solid's ultra fast dielectric response to a moving charge 
with unprecedented accuracy. This is supported by our 



quantitative prediction of wake-induced delays exceeding 
50 as that fall within the temporal resolution achievable 
with contemporary laser technology [lo| . 
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